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I. INTRODUCTION 

Brane-world models with a single extra dimension mm are decidedly a 5D phenomenological 
realization of Hofava-Witten supergravity solutions [3], if the moduli effects from compact extra 
dimensions can be ignored (for a review, see e. g. [Ij). The Hofava-Witten solution [3] can be 
thought of as being effectively 5D, with an extra dimension that can be large, when collated to 
the fundamental scale. They provide the basis for the well-known Randall-Sundrum (RS) brane- 
world models mm, that comprise the mirror symmetry and a brane with tension as well, which 
counter balances the leverage of the negative bulk cosmological constant on the brane, encompassing 
furthermore the branes self-gravity [Ij. In RS brane-world scenarios, our Universe is embedded in 
a 5D bulk of type AdSs [2]. The formalism to be used hereon employs a general metric for the 
brane-world, instead of the Minkowski metric in the standard RS model [2] . 

Brane-world black holes were comprehensively studied in Randall-Sundrum like brane-world 
cosmologies m\, where the dynamical equations on the brane are different from the general 
relativity ones. In fact, the brane-world framework presents terms that handle both the effects of 
the free gravitational field in the bulk and of the brane embedding in the bulk as well. The imprint of 
the nonlocal gravitational field in the bulk on the brane provides a splitting into anisotropic stress, 
flux, and nonlocal energy density, where this last determines the tidal acceleration out of the brane, 
possibly opposing the formation of singularities [3]. Unlike the nonlocal energy density and flux, the 
nonlocal anisotropic stress is not ascertained by any evolution equation on the brane. In particular, 
isotropy of the cosmic microwave background make the existence of the FRW background be under 
risk. Adiabatic density perturbations are furthermore coupled to perturbations in the bulk field, 
making an open system on the brane [7]. 

Bulk effects in the cosmological dynamics of brane-world scenarios have been studied [8] , also 
in the context of thick brane-worlds [M3], brane-worlds with variable tension [HHni and Kalb- 
Ramond fields [18], as well as in the context of more general brane-worlds and bulks |19H23j . 
Moreover, consequences of the gravitational collapse were proposed in the context of brane-world 
scenarios in, e. g., [23H3H]- In addition, dark matter was investigated already in [33] as a bulk 
effect on the brane. 

Black strings can be thought as being extended objects endowed with an event horizon, in low 
energy string theory m- The bulk metric near the brane and the black string warped horizon 
along the extra dimension are here reviewed, based on previous developments mm- Originally, a 
Schwarzschild black hole on the brane-world was shown to be a black string in a higher dimensional 
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spacetime, what leads to the usual astrophysical properties of black holes to be recovered in this 
scenario m- In this prototypical context, the Kretschmann curvature invariants diverge when the 
black string event horizon is approached along the axis of the black string. Several generaliza¬ 
tions provide attempts to preclude singularities both in the bulk and on the brane-world as well. 
When variable brane tension scenarios are taken into account, the brane tension can gauge the 
Kretschmann scalars involved. For instance, regular bulk solutions and black strings were obtained 
in Friedmann-Robertson-Walker brane-worlds under the Eotvos law [32], where the singularities 
related to the McVittie metric can be partially controlled as the cosmological time elapses. Indeed, 
for this type of metric the 5D physical soft singularities in the bulk are alleviated as time elapses, 
providing a regular 5D bulk solution, as the 5D Kretschmann invariants do not diverge. When 
other metrics are taken into account, for instance the Casadio-Fabbri-Mazzacurati metric one [33], 
black strings can be still emulated [31], however the related singularities in the bulk remain, re¬ 
gardless. In order to accomplish it, effective/perturbative approaches are usually employed, where 
the black string is made to evolve from the brane-world [IIEIIISS]. 

Brane solutions of static black hole exteriors with 5D corrections to the Schwarzschild metric 
have been found, for instance, in |36H38| . and furthermore in the context where the bulk singularities 
can be removed [39]. The (Schwarzchild) black string is unstable near the AdSs horizon, defining 
the so called Gregory-Laflamme instability iaiii]. This scenario might be drastically altered 
by the inhomogeneous dust and the dark radiation. In order to accomplish this effect, we use 
a procedure to calculate both the metric near the brane and the 5D black string horizon [15] . 
uniquely from a brane-world black hole metric and the associated Weyl tensor. Based on the 
knowledge of both the Sasaki-Shiromizu-Maeda effective field equations on the brane and upon the 
5D Einstein and Bianchi equations [H [T] 135] H2l 03] , both the bulk metric near the brane and in 
particular the black string warped horizon can be designed, by using a Taylor expansion along the 
extra dimension. Such procedure provides information about all the bulk metric components |15j . 
Indeed, the bulk spacetime may be either given, by solving the full 5D equations or alternatively 
obtained by evolving the brane-world black hole metric off the brane, what encompasses the imprint 
from the bulk via the Weyl tensor. Numerical methods have been employed to find black hole 
solutions in the context of black strings and fluid/gravity correspondence [H]. Similar methods 
involving expansions of the metric have been used in the context of black strings |44j . disposing 
the black string metric as the leading order solution in a Taylor expansion. 

The bulk shape of the black string horizon has been merely investigated in very particular 
cases mm, and latterly the standard black string was studied in the context of a brane-world 
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with variable tension m- Moreover, realistic models that take into account a post-Newtonian pa¬ 
rameter on the Casadio-Fabbri-Mazzacurati black string |45) , and the black string in a Friedmann- 
Robertson-Walker Eotvos brane-world |34j . also represent interesting applications. 

Recently regular black strings solutions associated to a dynamical brane-world have been ob¬ 
tained in the context of a variable brane tension [32] . The analysis of the 5D Kretschmann invariants 
makes us capable of attenuating the bulk physical singularities along some eras of the evolution 
of the Universe for the McVittie metric on an Eotvos fluid brane-world. This paper is devoted to 
encompass a framework with dark radiation and inhomogeneous dust. The 5D physical singulari¬ 
ties in the bulk are shown to be inherited from the 4D brane-world and no additional singularity 
appears in the bulk, for some range of parameters in our model. Nevertheless, the bulk physical 
soft singularities can be unexpectedly controlled in the bulk upon time evolution, what makes a 
regular bulk 5D solution in most ranges of the dark radiation parameter. 

This paper is organized as follows: in Section II the dark radiation dynamics on the brane is 
analyzed and reviewed. Starting with the Lemaitre-Tolman-Bondi (LTB) metric on the brane, 
the effective field eqnations for dark radiation on the brane are solved. The dynamical radiation 
model is shown to mimic a 4D cosmological constant on the brane. Both the black string solution 
and the bulk metric are obtained thereon. After obtaining the standard dark radiation model, 
a generalized framework is proposed. Both associated metrics are derived. In Section III the 
5D bulk metric near the brane and the generalized black string are derived and studied, and in 
Section IV the black string warped horizon in the context of inhomogeneous dust and generalized 
dark radiation is studied. Moreover, the black string physical singularities are analyzed from the 
Kretschmann invariants. The 5D physical singularities in the bulk reflect the 4D brane-world 
physical singularities. We analyze further the Kretschmann scalars generated by higher order 
derivatives of the Riemann tensor, and the respective physical soft singularities show that the bulk 
5D solution is regular, in some ranges of the dark radiation parameter. 

In order to fix the notation hereupon /i, = 0,1, 2, 3 and M, N = 0,1, 2, 3, 5, and let n be a 

time-like covector field normal to the brane and y the associated Gaussian coordinate. The brane 
metric components and the corresponding components of the bulk metric are related by 
gfiu + nf^riu = g^iv |1|- With these choices we can write 555 = 1 and 5^5 = 0, and thus the 5D bulk 
metric reads 

guN dx^ dx^ = g^u{x'^,y) dx^ dx'" + dy'^ , (1) 


where M, V = 0,1, 2, 3 effectively. 
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The initial paradigm concerning a perturbative method for obtaining the black string solution 
consisting in assuming the Schwarzschild form for the induced brane metric, on a RS brane-world. 
Subsequently a sheaf of such solutions are disposed into the extra dimension Eni: 


where ^ — \J^ denotes the curvature radius of the bulk AdSs, wherein the RS brane-world is 
embedded. Each space of constant y is a 4D Schwarzschild spacetime with a singularity along 
r = 0 for all y, the well-known (Schwarzschild) black string. 

As it is going to be clearer in Section III, the areal radius of the sheaf of such solutions along the 
extra dimension is called the black string warp horizon, that shall be precisely defined in Section 
III. 


II. DARK RADIATION DYNAMICS ON THE BRANE 

Henceforward some results concerning the dark radiation dynamics on the brane will be shortly 
revisited |46H48] . in order to briefly introduce the framework to get both the bulk metric near the 
brane and the black string encompassing the dark radiation parameter and the effective cosmolog¬ 
ical constant. New black strings solutions are here derived in the scenario provided in |46H48j on a 
background described by a Lemaitre-Tolman-Bondi (LTB) metric |49H51] . A solution for the black 
string that has as limit a tidal Reissner-Nordstrom black hole solution on the brane was obtained 
in [36], in the Randall-Sundrum scenario. In order to work with the effective Einstein equations on 
the brane, some conditions on the projected Weyl tensor are usually assumed, in order to provide 
a closed system. Besides, by taking into account a system of equations where a specific state equa¬ 
tion leads to a inhomogeneous density that has precisely the dark radiation form |52H54j (and its 
generalizations for thick branes OESI), and by solving the effective Einstein equations, the LTB 
metric can be derived. Subsequently both the bulk near the brane and the black string solution is 
obtained are obtained in this context. 

Erom the dynamics of a brane-world with matter fields present in the bulk m. the associated 
black string solution will be shown to present a generalized dark radiation form. The way to obtain 
the LTB metric is essentially different from that acquired in |36|. The inhomogeneous density is 
associated with conformal bulk fields, instead of being related to the electric part of the Weyl 
tensor. The black string solution is similarly obtained by a change in the coordinate system and 
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its final form generalizes the first case. Thereafter two new black string solutions will be presented 
and subsequently used in the construction of the horizon profile in the bulk. 


A. The Lemaitre-Tolman-Bondi Metric 


The Lemaitre-Tolman-Bondi (LTB) metrics are exact solutions of the Einstein equations that 
describe inhomogeneous spacetimes, having dust as the source. This type of models consider 
inhomogeneous generalizations of the Friedmann-Robertson-Walker (FRW) metrics, as the LTB 
metrics. An alternative approach for the LTB space-time is based in evolution equations of co¬ 
variant objects, as the density, expansion scalar, electric Weyl tensor, shear tensor and spatial 
curvature. The dynamics is reduced to scalar equations, and the FRW spacetime is achieved when 
two of these scalars associated with the shear tensor and electric Weyl tensor are zero. This formu¬ 
lation is based on a 1 -|- 3 covariant description [56] , which can be further applied to the LTB model 
|57l|58|. In general, the applications of these models involve black holes, galaxy clusters, super¬ 
clusters, cosmic voids, supernovas and redshift drift, for instance jSSj- Initially found by Lemaitre 
|49j , the LTB metric describes a spherically symmetric inhomogeneous fluid with anisotropic pres- 
surecosmological constant are present for instance in the Tolman model |50|. To derive the LTB 
solution, in comoving coordinates the general form for the line element is given by: 


= gfjLudx^dx'' = —dt^ -|- A^{r,t)dr‘^ + E?{r,t)dQ.2, 


(3) 


where the 2 -sphere area element is denoted by dQ- 2 , the energy-momentum tensor written as Tj^ = 
diag{—p, —pA, —PA, —Pa), where p denotes the energy density, A 4 denotes the brane cosmological 
constant with associated energy density pA = «: 4 ^A 4 , and K 4 is the 4D gravitational coupling 
constant. The Einstein equations, for each one of the space diagonal components, are given by the 
following expressions: 


1 

R 



1 + {dtKf - 


d^R dtR dtA drR dr A dfA 

~ir ^ ^ ^ ~T~ 


drR 


Rj^2^tdrR 

d^ 

A 


drR 


A4, 

A4, 

A4. 


(4) 

(5) 

( 6 ) 


The function A{r,t) = g{r)drR satisfies Eq.([^. By setting g{r) = (1 -|- /(r)) the usual form 
of the LTB metric is hence obtained: 



-dt^ + dr"^ + R^dn2, 


1 + / 


(7) 
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where /(r) > —1. The function / can be interpreted as the energy density shell f{r) = 2E{r). The 
function g{r) is a geometric factor such that when g{r) = 1 the spatial sections are flat. Eqs.Q 
and (j^ are not independent, leading to the expression 

{dtRf = f + ^ + ^R^, ( 8 ) 


where M = M{r) is an arbitrary function of integration that gives the gravitational mass within 
each comoving shell of coordinate radius r. By definition of mass in one can write 2dm/dr = 
k\pAR?, which implies M = J drm\/\ + fdr. The non-relativistic limit / <C 1 yields M ^ m, and 
Eq.([^ reads ^ {dtR)^ ~ ^ ~ 2 ' term is interpreted as kinetic energy, the second stands 

for the Newtonian potential term and / is twice the energy of the system, when A 4 = 0. Hence 
M is the relativistic generalization of the Newtonian mass. When / is negligible, namely in the 
non-relativistic limit, the spatial sections are flat when g = 1. The function g provides the energy 
in each spatial section, and thus carries the information of curvature for each section. 

Finally, Eq. Q implies that 


t - tN{r) = 


cR 


where tpf is known as the “bang time”. Eq. 
classes. When A 4 = 0 it reads: 


/ + ^ + iR, ( 9 ) 

(fsl) can be used to classify the LTB models into three 


-l </<0 


/ = 0 


/>0 


elliptic, 
parabolic, 
hyperbolic. 


When A 4 7 ^ 0 the potential V{R) = ^ + ^R? leads to a different classification, depending on 
the sign of A 4 . 


B. LTB Solution on the Brane 


In this Subsection the LTB solution associated with dark radiation on the brane is reviewed, 
starting with the projected Einstein equations on the brane in vacuum. Unlike the Reissner- 
Nordstrom black hole, this new solution has a specihc dark radiation tidal charge Q. The 4D and 
5D coupling constants are related by k\ = The field equations in 5D Einstein theory lead 

to the projected equations [l 2 l [60] 


G 


flU — 








( 10 ) 
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where is a term quadratic in the energy-momentum and provides high-energy corrections 
arising from the extrinsic curvature of the brane, what increases the pressure and effective density 
of collapsing matter. The term is the projection of the bulk Weyl tensor and provide Kaluza- 
Klein corrections originated from 5D graviton stresses [28], as the massive modes for the graviton 
in the linearized regime. For observers on the brane such stresses are nonlocal, in the sense 
that they are local density inhomogeneities on the brane generate Weyl curvature in the bulk 
backreacting nonlocally on the brane mizi HI iniiHa [51]. Therefore the vacuum equations are 
where the Weyl projected tensor can be identified with a trace-free energy- 
provided by [HI] 

^ r / T \ 1 

( 11 ) 


momentum as ~ 


f — 


U ( ^ h^i/ I + 


where U is the effective energy, P^i, is an anisotropic stress tensor, is the effective energy flux, 
is a 4D velocity vector satisfying = —1, and is such that = 0, being thus possible 

to write + v^j^Vy. A non-static spherically symmetric brane-world with P^i, ^ 0 can be 

described by the line element 


'^'5(4) = t)dt^ + B^{r, t)dr^ -|- t)dEl 2 


( 12 ) 


with Qfj_ = 0. The anisotropic stress tensor can be represented by P^u = P (r^rj/ -|- 
where P = P{r,t) is a scalar field and is the unit radial vector. With these assumptions the 
electric part of the Weyl tensor yields 


= { — ] diag{p, -pr, -pr, -pr), 




(13) 


with p = U, 3pr = U + 2P and 3pT = U — P. 

By assuming the brane field equation = 0 |1] and by considering the state equation 

, implying that 


p = -Pr, it yields dtU + = 0 = + A^U 


U=(^A ^ 

where the constant Q is the dark radiation tidal charge 


(14) 

As p = U, the energy density in 


this case is related to the inhomogeneous density. Thus the 4D Einstein equations (10) read 

G = Ai^g^n 2ry P h^y). 


(15) 


It follows, by solving the above equations, that the component Gtr is obtained by the expression 
dtBdrR — BdtdrR = 0, and therefore the function 


B = R-^drR 


( 16 ) 
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satisfies this relation with H = H{r). Considering thus such expression for B in the trace equation 

it is possible to write the expression as follows 


-G/ + GI + 2G!! 




(17) 


which is similar to Eq.([^. Integrating Eq. ( |17[ ), it reads 

±* + T(r) = y'(/-;| + l^) dR. 

which is analogous to Eq. ([^. It is thus possible to write (12) in the LTB form given by ([^. 
Making the transformation of the LTB coordinates (t, r) to curvature coordinates (T, R) as 


(18) 


T = t + 


the following 4D metric is finally obtained: 


-R2_Q 


dR, 


Q A4 


dsf^) = - I 1 + ^ - -^R" ]dT^+{l + ^- -^R 




Q A4 




-1 


dR + R^dQ,2 ■ 


(19) 


( 20 ) 


This metric is known as the inhomogeneous static exterior of a collapsing sphere of homogeneous 
dark radiation j28l I62| . Note that when A 4 = 0, the solution (20) is formally analogous to the 
Reissner-Nordstrom solution, when one identifies the electric charge to the dark radiation tidal 
charge. 

In what follows this solution will be generalized, by considering a generalized dark radiation 
term with dark radiation charge Qn where 7 / is a parameter characterizing the model of the dark 
radiation m- The dynamics of a spherically symmetric brane-world is also analyzed, when the bulk 
a) carries matter fields; and b) when its warp factor characterizes a global conformal transformation 
consistent with Z 2 symmetry. Einally, it is possible to study the bulk metric and the black string 
solutions with a term analogous to the black hole solution with cosmological constant on the brane. 
In this framework, the energy-momentum tensor encompasses conformal bulk matter fields, whose 
dynamics provide a specific state equation [IGIHH] . 

Consider now a general conformal spherically symmetric metric ds\ consistent with Z 2 symmetry 
along the extra dimension on the brane 


dsl = + e^^dr'^ -|- R?dVL2 + dz^) , 


( 21 ) 


where z stands for the conformal extra dimensional coordinate, and A, B, R, and 11 are general 
functions of the coordinates {t,r,z). H denotes the conformal factor. The Einstein field equations 
are given by 

Gmn = —kI —Tmn + ^d{z — zo)gMN + , (22) 
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where guN denotes the components of the is the induced metric and T^u stands for the components 
of the energy-momentum tensor representing the bulk fields. The brane is localized at 2 ; = fo¬ 
under the conformal transformation , the energy-momentum tensor is assumed, 

as usual, to have weight s = —4 [l6]. The conformal Einstein tensor is given by 


Gmn = Gmn + '^MNig) ■ 


By using the expression g^ =14 ^g^, Eq.(22) hence leads to the following equations: 


G M _ 

N — ^ 5 -^ N -> 

= - zo)g^ + ] • 


1 M 


M 


(23) 


(24) 

(25) 


Eq.(24) evinces that the 5D Einstein tensor is related solely to the presence of fields in the bulk. 


and is independent both on the brane tension A and upon the bulk cosmological constant A5. 
Moreover, Eq. (|25l) provides the dynamics of the conformal factor. The divergence condition leads 


to ^ — Td^^] = 0, and from the Bianchi identity Eqs. (24) and (25) it 

implies that 


either 


or 


VmT% = 0 
3T^dM^ = TOn^ . 


(26) 

(27) 


Hence Eqs. (26) and (27) require that 2T/ = Now, by considering the energy-momentum 
tensor 


= diag{-p,pr,PT,PT,Pz) , 


(28) 


the state equation p — Pr — 2pT + 2pz = 0 holds. It implies that VzT/ = 0 and subsequently that 
dzPz = 0 Thus p, Pr and pT must be independent of 2 . 

The system of inhomogeneous dark radiation and an effective cosmological constant is defined 
by a conformal bulk matter with the following equations of state 


A , A , , 

p = -Pr = PDR + and PT + gp = ^ {11 - 1) J 


(29) 


"'5 "'S 

where g characterizes the dark radiation model and A is a bulk quantity and mimics a 4D cosmo¬ 
logical constant on the brane. The components of the Einstein tensor can be thus evinced: 


g; = g/ = -k2p^r-a, 

G/ = —k^pdr {1 + g) — 2K . 


Gg^ = G^ = -4gpDR-A, 


(30) 

(31) 
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By taking the divergence of Tj^, it reads 


dtPDR + {1 - rj) = 0 = drPDR + ‘^PDR-^ (1 “ > 


(32) 


what consequently leads to 


where here Qr 


{ 14) when p = 


Pdr = (33) 

= const is interpreted as a generalized dark radiation tidal charge and behaves like 
— 1. Regarding the following energy conditions |63) 


p>0 

and 

P + Pi >0 

P + '^Pi > 0 

and 

P + Pi >0 


2=1 


P > \Pi\ , 


(34) 

(35) 

(36) 


respectively known as weak, strong and dominant energy conditions, for /? = 0 (by taking the weak 
energy condition) the equality poR = holds, such that A = < 0 for > 0. 

One can realize that p = —pr regards the weak condition and p + pT = {I — p)pdr holds if 
r/ < 1, where here /? > 0 satisfies the first condition. By a similar analysis in 5D we see that p <0 
(weak), and p < 0 (4D and strong), p < —1/3 (5D and strong), \p\ < 1 (4D and dominant) and 
—2 < p < 0 (5D and dominant). If < 0, all the energy conditions are violated. Taking the 
trace equation and then making two integrations, it yields 


idtRf = :^R^^ + ^R^ + f, for ry/-1/2, (37) 

2p + 1 3 

what generalizes Eq. ([iTl), being identical to it when p = —1. By performing one more integration 
it reads 


± t + r 


Qr] 

2p + l 






(38) 


where refers to expansion or collapse and r corresponds to the evaluation of the function at 
t = 0. The condition R{0, r) = r is taken on the hypersurface t = 0. The radial equation leads to 
77 = VT+TIIZ!, and thus the metric has the 4D LTB form. Hence the 5D conformal line element 
reads 


+ {drR/H)^ dr"^ + R^d^2 + dz^'^ , (39) 

where ^rs is the Randall-Sundrum warp factor [2j. From the dynamical dark radiation models the 
marginal bound / = 0 corresponds actually to static solutions. Finally, the transformation from 
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LTB coordinates (t, r) to the curvature coordinates (T, R) is given thus by 


T = t + 


Qv U2ri I A p2 
2 r)+l-^ “'“3'^ 


i?2?? + A P2 _ I 
2'q+l^ ^ 3 -^ 




(40) 


The black string solution is therefore obtained: 


'^'5(4) = - I 1 - 


2r/ + 1 3 


- ^—R2n _ ^r2 

2r] + l 3 


-1 


dR‘^ + R^dn2. 


(41) 


The models that describe the inhomogeneous static exterior of a collapsing sphere of homo¬ 
geneous standard dark radiation require the value rj = —1 |28[ 1461162j . In general the exterior 
spacetime is not static in the brane-world scenario [28], however the collapse of a homogeneous 
Kaluza-Klein energy density is static, and can be identified to the dark radiation. Hence in the 
case r] = —1 the model of homogeneous dark radiation is recovered. It is worth to emphasize that 
the exterior is static solely when the system has tidal charge and cosmological constant, where the 
physical mass equals zero. When A = 0, the zero mass limit of the tidal Reissner-Nordstrom black 
hole is obtained [36]. The event horizon is determined by the solutions of the following equation: 


1 - _ ^j^2 ^ Q (- 42 ) 

2ri + l 3 ^ ^ 

on the brane. Thus for A = 0 it implies that Rf^ = {2r]-\-l)/Qn. When A 7 ^ 0, the exact location of 
the horizons can not be obtained, with exceptions for the values rj = —1 and rj = 1/2. Respectively, 
the horizons are given by: 


„(±) ^ I 

I 2A 


d(±) _ ^*31/2 

- 4 ^ 


1 ± 1 + 


4Q_iAy/^' 


- 1 ± 1 + 


16A 


1/2 


9Q2 


1/2 

> , for 7 / = —1, 

, for r? = 1/2 . 


1 / 2 . 


(43) 

(44) 


The conditions A < 0 and Q 1/2 > 0 imply two horizons, an inner and an outer one. For two 
specific values of Q 1/2 the singularity is naked, thus violating the cosmic censorship hypothesis 
|64j . Yet, if A > 0 and Q 1/2 ^ 0, then there is a single horizon R^^'^ ■ 


III. BULK METRIC AND THE BLACK STRING 

In this Section, the bulk metric near the brane as well as the black string associated to a black 
hole on a brane-world are briefly introduced Eq.([^ represents the black string metric. 

By denoting the components of the 5D Riemann tensor, as the 5D Kretschmann invariant 
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*'^^’^Q/3p(T g4|p|/£ jg unbounded as y —>• oo [H [3T], thus the Schwarzschild 

solution is not a good candidate neither for a brane-world black hole nor for trying to remove at 
least some of the bulk singularities. Hence a well-established perturbative method is employed 
to find both the bulk metric near the brane and, in particular, the black string warped horizon 
along the extra dimension. In what follows such framework is revisited. This shall be further 
accomplished in Section IV in the context of inhomogeneous dust and generalized dark radiation 
on the brane-world and inhomogeneous dark radiation in the bulk as well. There we shall evince 
that in such scenario the bulk (and the black string) can be regular for certain ranges of the dark 
radiation parameters. 

In brane-worlds with Z 2 symmetry, the junction conditions imply that the extrinsic curvature 
of the brane is given by |1] 


K 


2 r 


T^i/ ^ (T A) 


(45) 


The trace-free and symmetric components of the bulk Weyl tensor C^yap are respectively given by 
^pLvp = dp Qy Crapis and Spy = Cpy^p n”" n^, where n“ denote components of a vector field out 
of the braneworld. Hereupon we denote by TZpy^p = ^^^'R-pyapix’^,0) the components of the 5D 
Riemann tensor computed on the brane. 

The effective field equations encompass the 5D Bianchi and Einstein equations [U [23l H2] , which 
have been generalized in a variable brane tension framework as well [IllliailllESlES]. The effective 
field equations are provided by the following system of equations: 


£ Kpy = -^A^gpy - Spy + KppK^y (46) 

£Sp. = ^A5{Kpy-Kgpy)+V^Br(py)+'Rppy„KP^+Kp^pK^}yKP^+3KypSy)p-KSpy (47) 

£ Bpyp = Kp^^Bpya - 2'SIypSy^p - 2Bp^[pKyf (48) 

£'T^pvpa = ~‘2£i,pvT[p£^cr\ ~'^p^[p(j]u- 


They are worked out the boundary condition Bpyp = on the brane [3j. These expressions 

can be employed to compute a Taylor expansion of the brane-world metric along the extra dimen¬ 
sion, what perturbatively defines the bulk metric near the brane and the black string horizon. The 
set of equations obtained from the 5D Einstein and Bianchi equations by Shiromizu, Sasaki and 
Maeda mils] provides the bulk metric near the brane, given by the standard Taylor expansion 
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along the extra dimension y: 


y) = giMuix"", 0 ) + {£ gi,u{x°‘, y)) \y=o \y\ + {£ {£ gf,u{x°‘,y))) \y=o ^ 
H + {£^{gyuix°', y))\y=o^-^ + • • • 


(50) 


As comprehensively studied in m, the term in the first order \y\ above is obtained from the 
definition of the extrinsic curvature = \£gy,u and from the junction condition (45) as well, by 


denoting g^i, = gy^{x°‘,0) hereon. The term in is proportional to £Kyy, at Eq.(46). Moreover, 
the term Ky_aK°‘v in Eq.(46) can be computed by (45). Besides, the coefficient term of |yp in 


Eq. (50) can be expressed as 




2 £ {£K^,) = £ -££y,- -^£g, 


6 


I ) 


(51) 


in virtue of Eq. (46). As the Lie derivatives terms in the right hand side of this last expression are 


respectively given by Eq.(46), due to Eq.(47), and by the definition of the extrinsic curvature, one 
hence arrives at the expression for |y|^ in Eq.(50). Einally, the term in y^ is acquired by a similar 
reasoning. 

One thus finds the above expansion up to the fourth order is given by m 


g^iv{x ,y) — g^i/ Kg 
+ ■ 


+ 


2 (A T)gy^u + 


\y\ 


-4 -(A - T)T^, + - - Ag - -4{\ - ry 


2Kyf,KiKy - + ^Agy^^iL + - K£, 


\u) ■ -pyiJ-i'j ' g 


2 ! 




l?/P 

3! 


As 

T 


+ 

t 

+ ( 


(n-^ + gy, + - As) Ky^Kl + (7^ - Ag + 2Ky£y, 

+ TZyapg + K Ky, - ^ Ag7^^, + 2KyfiK%KP^Ky 


+ ]-KKy£,^ - lK^^Ky^n,^a0 + Spa ( \KKy - 
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--Ky^sixy - 


2 


4! 


(52) 


where H = Hy and = Hp^jHP^, for any tensor H of rank two. 


The black string warped horizon [66] can be studied from gg 0 (x'^,y) in Eq. (52). In fact, 
regarding a spherically symmetric 4D metric 


'^' 5 ^ 4 ) = —hi{r) dy + h2{r) d4 + 4 d£l2 , 


( 53 ) 
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modelling a brane black hole, the usual 4D areal radial coordinate r is related to the 5D metric 
in Eq.(52) by the expression = g 0 g{x‘^,O) |15[ I35| . The black string has a warped horizon on 
the brane that has radius r = gg 0 {x°‘, 0), where r denotes the coordinate singularity which is 

solution of grr{ic) = 0 (See Eq. (II.6) of Ref. [67]). The radius of the black string warped horizon 
is hence r(y) = \/g 00 {x‘^, y). The term g 00 {x^, y) in Eq. (52) for g = u = 9 corresponds to the bulk 
(squared) areal radius, and includes both the black string horizon for r = r, and in particular the 
brane black hole horizon for y = 0 . 


IV. BLACK STRINGS AND DARK DUST 


In this Section we investigate the black string related to the induced black hole on the brane 


given by Eq. (41). The black string warped horizon is well known to be the component g 00 {x°‘,y) 


in Eq.(52), evaluated at event horizons in (44). At first let us calculate such component for an 
arbitrary x“. In order to accomplish it, the 4D energy-momentum tensor on the brane is given by 


the brane-world components in Eq.(28), where 


P = -Pr = ^ ^ 

(54) 



ni ni 

(55) 


and the expressions for the energy-momentum, the projected Weyl tensor components, the extrinsic 
curvature, and the Riemann tensor respectively given respectively by ( Ala||A4c) in Appendix are 


thus obtained. Hence the component 00 of the metric, corresponding to the black string horizon 
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along the bulk (52), can be written as: 


9ed{x‘^,y) ^ (4 - r/) + 3A] |y| 


^ (1 — r/) + - (Ak| + 2 ^(?7 + 2 ) + 3A]) — A 5 

“2 ^ + A) [2Ak| + QrjR^'^ ^ (r/ + 8 ) + 9A] | — 


+R 


1^q^r2v-2 (^ _ 4) _ _ 3 ^] 


-^As [QrjR^^-^ (r? - 4) - AkI - 3A] - ^R^v-2 (4 _ [2AKi + Q^i? 2,-2 ^ 2 )] 


+Q,R^^-^ (r? - 1 ) [^R^^ ( 7 ? - 1 ) + 1 +ga 2 , 4 _ 


Qr 


A 


-1 


- __j^2 j exp(-aT) 

2r? + l 3 ' ^ 

^ r2 

■3^ 


( 2 ^ - 2) (77 + 2) (27? - 3) (^1 - 
+ 3 ^As - Ak^ - A] - 1 f+ A] 

+^R^v-2 (1 _ [q^r2v-2 (^ _ 4) - AkI - 3A] 

216 ^ “ 4) — Ak| — 3A] ^ (977 — 5) + 3Ak| — 3A] | LI |_ ... 


(56) 


The last expression does not explicitly include the terms of order 7 /^/ 4 ! as they are extensive and 
awkward, although we shall consider such terms in our subsequent analysis. 

Our first analysis consists in regarding the variable T = T(r,t), which appears solely as the 
term exp(—aT) at the order | 7 /p/ 3 ! in Eq. (56). The graphic below in Fig. [^shows the dependence 
of the component gee{x°‘,y) on T = T{r,t). Hence, near the brane, our results can be taken to be 


<7 = 1/2 



FIG. 1. Black string warped horizon gge(R,y) along the extra dimension for A = — 1, for a = 1. 


independent of the coordinate T. Below in Figs. EE and 1 ^ we depict the black string warped 
horizon ggg{R,y) as a function of the parameter r/, for the possible values A = 0, ±1: 
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A= 1 



FIG. 2. Black string warped horizon area 
gee {R, y) along the extra dimension for 
A = 1 and the parameter g variable. 


A = 0 



FIG. 3. Black string warped horizon area 
gee {R, y) along the extra dimension for 
A = 0 and the parameter g variable. 



A = -l 


FIG. 4. Black string warped horizon area gee{R,y) along the extra dimension for A = — 1 and the parameter g 
variable. 


Fig. ^ illustrates the case of a positive cosmological constant A = 1, where the dark radiation 
parameter ij plays an important role on the black string warped horizon profile. The bigger the 
value of rj, the faster the warped horizon goes to zero along the extra dimension. In the range 
1-8 ^ ^ 8.7 the black string warped horizon increases slower along the extra dimension, however 

it still goes to zero for higher values of the extra dimension. Figs. and show, respectively 
for A = 0 and A = — 1 , quite different prohles, as the black string warped horizon is shown 
monotonically to increase along the extra dimension, as rj increases. 

Now, the black string warped horizon is analyzed for an inhomogeneous static exterior of a 
collapsing sphere with homogeneous standard dark radiation, as a function of Qr^, as we have seen 
that when A 7 ^: 0 in general it is not possible to obtain the exact location of the horizons by 


Eq.(42). In fact the two single exceptions are the models corresponding to r/ = — 1 and rj = 1/2, 
given by Eqs.(44). For fixed value of the dark radiation parameter. Figs. and show that the 
generalized dark radiation tidal charge Qn also influences the black string warped horizon along 
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ri = -l 



FIG. 5. Black string warped hori¬ 

zon area gBe{R,y) along the extra dimen¬ 
sion, as a function of the generalized dark 
radiation tidal charge parameter Qn, for 
rj = -1. 


77 = 1/2 



FIG. 6. Black string warped horizon 
area gge{R,y) along the extra dimension, 
as a function of for g = 1/2. 


the extra dimension. When r/ = —1, Fig. [^illustrates that the horizon goes to zero for some point 
along the extra dimension for Q-i > 5.1, while Fig. [^evinces that the horizon goes to zero for 
some point along the extra dimension for Q 1/2 ^ 1-2. 


Besides the black string warped horizon, Eq.(52) provides more generally the bulk metric near 
the brane. In order to check whether the bulk is regular, we aim to study the 4D and the 5D 
Kretschmann invariants, related to the black hole on the brane and to the black string in the bulk. 


Our analysis of the bulk physical singularities is independent on the perturbative method (52), as 
the curvature invariants are independent of it. 


The Gauss equation which relates the 4D and 5D Riemann curvature tensors according to 


and the 4D and 5D Kretschmann invariants are given by 

k(4) := 

^( 5 ) _ 


(57) 

(58) 

(59) 


For the generalized dark radiation model they are related by: 
Ki^) = ^(4) _ ^2 (ry _ 4) - - 3A] 


_ 1 ) 




+ — 4 

18 I 


-AkI -|- ^ (77 — 1) 


n 2 




(60) 
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The physical singularities are provided by i? = 0 and by the solutions Rh of the Eq. (42) for the 


values rj = —1 and rj = 1/2 given by Eqs.(44). Neither new singularities are introduced nor the 


existing ones (i? = 0 and R = R^ in Eqs.(43) and (44) are removed from the bulk. 

We are going to show in what follows, by the analysis of the physical soft singularities, that the 
bulk can be regular. Indeed, the 4D invariant 




(61) 


(here denotes the covariant derivative on the brane) is very soft, since it takes invariants 
involving at least two derivatives of the curvature to detect it. Consequently, the 5D version of the 
invariant ^ reads 

(5)^ = (62) 


where Da denotes the 5D covariant derivative. It is worth to point that a, b are effectively 4D 
spacetime indexes, as the 5D covariant derivative can be realized as Da = and Da = V 5 , when 
the extra dimension y is taken into account. It implies that the 5D Kretschmann invariant is 
given by [l5] 


-^{VlKrpK^a){m^R^^n + 2{VlK,ya\K^\p^){m'^K^^K^n • (63) 


Based on the values of the metric (41) and the extrinsic curvature components given in (45), the 5D 
Kretschmann invariants for the bulk can calculated. Due to the awkwardness of the expression 
for these invariants, we opt to analyze our results by the graphics in Eigs. 8, 10, 12, and 14. Eig. 
[^has to be compared to Fig. that describes the 5D invariant In Fig. since r/ ~ 0 then 
the 4D invariant ^ diverges, independently of the value for R. On the other hand, Fig. [^depicts 
that for 7 / ~ 0, the 5D Kretschmann scale goes to infinity for most values of i?, but when 
0.94 ^ R ^ 1.06 the 5D Kretschmann invariant does not diverge. Hence for this range of 
physical soft singularities present on the brane the bulk is regular. Moreover, Fig. shows that 
^ — >■ +00 for i? —)• 0, however Fig. [^evinces that for 0 < r/ < 0.5 the limit ^ —)> +00 on the brane 
alters to — )> —00, for values of ii ~ 0. Still, for —1 < < —0.5 the 5D Kretschmann invariant 

does not diverge for 0.08 ^ R ^ 0.13, and again this range of physical soft singularities on the 
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FIG. 7. 4D Kretschmann invariant ^ 
for Qn = 1, for A = 0. 


FIG. 8. 5D Kretschmann invariant for 
Q,, = 1, for A = 0. 


brane correspond to a regular bulk. Notwithstanding, now there is a range —1 < < —0.5 and 

0.08 — 0.13 ^ R 0.47 — 0.49 that makes the 5D Kretschmann invariant to diverge. 

Now the case A = 1 is regarded. Fig. regarding the 4D invariant must be compared 




FIG. 9. 4D Kretschmann invariant ^ 
for Qn = 1, for A = 1. 


FIG. 10. 5D Kretschmann invariant 
for Qn = 1, for A = 1. 


respectively to Fig. [lol describing the 5D invariant A similar pattern is realized in this case, 
where now Fig. shows that ^ —)> +oo for i? —)• 0. We regard the A = — 1 case in Figs. 11 and 12 


and the conclusions are analogous to both previous cases. Finally, the values of the 4D and the 5D 
Kretschmann invariants are provided for = 10 and A = 0. The respective graphics for A = ±1 
are quite similar. 


Fig. 13 illustrates that whatever the value for r] in the range [—1,1/2], the 5D curvature 
invariant is finite for R > 0.71, and the whole bulk is regular, independent of the value of the 


extra dimension y. Nevertheless, Fig. 14 also evinces the same pattern of the previous cases. For 
0 < 7/ < 0.5 the limit ^ —?■ +00 on the brane is modified into the bulk to —)■ — 00 , for values of 
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FIG. 11. 4D Kretschmann invariant 
^ for = 1, for A = —1. 


FIG. 12. 5D Kretschmann invariant 
(5)^ for = 1, for A = -1. 




e, = 10 

0 


FIG. 13. 4D Kretschmann invariant ^ 
for Qn = 10, for A = 0. 


FIG. 14. 5D Kretschmann invariant 
for Qr, — 10, for A = 0. 


R near R = 0. Still, again for —1 ^ rj < —0.5 the 5D Kretschmann invariant does not diverge 
for 0.08 < ii < 0.13, and again this range of physical soft singularities, present on the brane, are 
banished in the bulk, as the Kretschmann invariants remain finite from a starting hnite value for 
R. 


V. CONCLUDING REMARKS 

The bulk metric near the brane was obtained and in particular the black string warped horizon, 
for inhomogeneous dust and generalized dark radiation on the brane-world and inhomogeneous 
dark radiation in the bulk are regarded. The standard dark radiation [ISIESHSI] is a particular 
case analyzed, where the dark radiation parameter 77 = —1 and the mimicked cosmological constant 
A equals zero, provided by the 5D pressure also equal to zero. 


By analyzing the 4D and 5D standard Kretschmann invariants, respectively defined by Eqs. (58) 
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and (59), the Gauss equation is shown to imply that the bulk associated to brane-world models with 


inhomogeneous dust and generalized dark radiation inherits the brane-world physical singularities 
at R = 0 and at i? = Rh is the solution of Eq. ( |42[ ). Although the black string warped horizon 
is obtained as a particular case of the 5D bulk metric near the brane via an effective method, the 
analysis of the 5D bulk singularities relies on an exact method provided by the Gauss equation 


(57) 


The 4D and 5D Kretschmann invariants, respectively given by (61) and (62), show that although 


some physical soft singularities exist on the brane, the 5D bulk can be regular, for some parameters 
of the dark radiation parameter, for various values of the generalized dark radiation tidal charge 
parameter Qrj analyzed. Figs. 7-14 show these features for some values of Qn and all possible 
values of the mimicked cosmological constant on the brane. In fact Figs. 7, 9, 11, and 13 depict 
the 4D Kretschmann invariants, while Figs. 8, 10, 12, and 14 illustrate the 5D ones, for all values 
of the mimicked 4D cosmological constant on the brane A = 0, ±1. Such figures evince that the 
bulk is more regular in what concerns the respective soft singularities, for certain ranges of the 
generalized dark radiation tidal charge parameter analyzed. 
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Appendix A: Energy-momentum, Weyl tensor and extrinsic curvature components 


The components of the energy-momentum and the Weyl tensor are given respectively by 






T0e = -2 + A) , = -2 + A) sin^ 0 (Ala) 

'^5 '^5 

R^’’ - + A), (Alb) 




2i( + l 




A 


-1 




(Ale) 
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See ={I - r]) 

0 


^iDip — ^QQ sin B , 


A 


-9i^jpv _ tPip 

2r/ +1 3 


Srr = ^{1-7^)R^^-^1- 


Qr, 


A 


- "^R^^ - -R^ 

277 + 1 3 


-1 


For the extrinsic curvature, by computing (45) we obtain 
Kee = Y (r? - 4) - A^i - 3A] , 


SSifif — R-99 sin 0^ 


Kj'j' = — ( 1 — 


Qr 


277 + 1 ' 


-R^^ - -R 


-\n\ + QrfR^^ ^ “ 1 ) 


Krr ^ “ 3 1 ^ “ 


_3^R2e _ 

277 + 1 3 


-1 r 


-Ak| + Qr^R^^ ^ (77 — 1 ) 


Finally, the necessary components of the Riemann tensor are given by 

Qn A 


Rerer = 


-_^2 _ ^ 

277 + 1 3 


RdRBR “ “ 

Re^e^p = R^ sin^ 


Qr. 


A 


-^R^r, _ _^2 

277 + 1 3 


-1 


Q]9j^J12v + 

277 + 1 3 

A 2 

277 + 1 3 


Qr 


A 


-^R 2 ^ + -R 2 

277 + 1 3 


(A 2 a) 

(A 2 b) 

(A2c) 

(A3a) 

(A3b) 

(A3c) 

(A4a) 

(A4b) 

(A4c) 
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